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SECTION  1 
INTRODUCTION 

1 . 1  PLASHOIDS 

Plasmoids  can  be  described  as  a  group  of  electrons  and  ions  that 
propagate  as  a  unit  through  vacuum.  In  their  present  conception  they  consist  of 
an  expanding  halo  region,  carrying  return  current  and  bearing  a  slight,  charge 
imbalance,  and  an  equilibrium  core  region  that  propagates  like  a  pinched  charged 
particle  beam.  This  model  was  first  suggested  by  T.  Lockner  of  Sandia.  These 
parts  are  illustrated  in  Figure  1.  Since  they  are  a  type  of  charged  particle 
beam,  they  can  be  generated  with  higher  efficiency  and  deposit  their  energy  over 
a  deeper  distance  in  matter  than  laser  beams.  This  makes  plasmoids  attractive 
for  defense  appl icat ions ;  however,  the  structure,  formation,  and  propagation  of 
plasmoids  in  a  space  environment  are  not  known  at  present.  Theory  and 
one-dimensional  simulations  aimed  at  finding  the  core  structure  of  plasmoids  by 
using  variational  principles  have  been  successful.  This  report  documents  the 
progress  made  in  this  first  phase  of  research.  A  proposed  second  phase  of 
research  will  use  similar  techniques  to  model  the  three-dimensional  formation 
and  propagation  of  piasmoids  in  a  space  environment. 

1.2  PROGRESS  IN  PHASE -ONE  RESEARCH 

The  first  phase  of  research  was  aimed  at  understanding  the  core 
equilibrium  of  a  plasmoid  as  being  the  result  of  a  relaxation  process.  In  this 
model  the  plasmoid  was  a  minimum  energy  structure,  and  internal  structure  could 
be  found  by  variational  calculus.  This  model  was  suggested  by  the  success  of 
such  methods  in  magnetic  fusion  work  (Reference  1). 

Progress  has  been  made  in  the  theory  of  plasmoid  core  equilibria  of 
three  types  (using  cylindrical  coordinates  with  z  being  the  equilibrium  axis): 

(1)  nonrotational  equilibria, 

(2)  rotational  equilibria  with  no  axial  magnetization 
(Vq^O  3y=0) ,  and 

(3)  rotational,  axially  magnetized  equilibria  (Vy*0  0^*0 
leading  to  a  Bz  field). 
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The  variational  method,  called  the  method  of  Lagrange  multipliers 
(Reference  2),  was  applied  originally  in  one  Lorentz  frame.  This  frame  of 
reference,  called  the  z-pinch  frame,  is  a  frame  where  the  core  equilibrium  is 
charqe  neutral  and  electrons  and  ions  move  past  each  other,  so  their  currents 
add.  The  plasmoid  core  equilibrium  can,  therefore,  be  understood  as  a  Lorentz 
transformed  z  pinch.  However,  the  variational  principal  originally  was 
formulated  in  the  z-pinch  frame  only  and  assumed  similar  density  profiles  for 
electrons  and  ions.  A  generalized  variational  principle,  valid  in  all  Lorentz 
frames  and  requiring  similar  profiles  for  both  electrons  and  ions  in  all  frames, 
has  now  been  formulated.  The  principal  of  Lorentz  invariance  strongly  indicated 
the  form  of  the  variational  problem  and  the  Lorentz  invariant  variational 
principal  was  found  to  predict  the  three  types  of  equilibria  mentioned  earlier. 

Simulations  of  the  three  types  of  equilibria  discussed  were  performed 
assuming  rotational  and  axial  symmetry,  that  is,  effectively  one-dimensional. 
The  simulations  consisted  of  equilibria  initialized  in  steady-state  and 
nonequilibrium  configurations  that  were  allowed  to  relax  to  equilibrium.  The 
simulations  demonstrated  that  the  equilibria  predicted  by  variational  theory  did 
exist  and  were  unique  end  states  of  the  relaxation  process.  Simulations  also 
demonstrated  that  relaxation  to  final  profiles  was  slower  than  the  expected  ion 
transit  time  timescale,  and  instead  proceeds  according  to  a  Bohm  diffusion 
timescale. 

1.3  PLASMOID  EQUILIBRIA  AMD  Z-PINCH  RAM- 

Before  proceeding  to  the  body  of  this  report  it  is  usefui  to  briefly 
consider  the  equations  of  plasmoid  equilibria.  A  simple  method  to  generate 
plasmoid  equilibria,  which  also  aids  in  understanding  them,  is  to  use  the  two 
dimensional  Virial  Theorem  (Reference  3).  This  theorem  gives  necessary 

conditions  for  the  existence  of  an  infinitely  long  rotationally  symmetric, 
pinched  equilibrium.  This  condition  is: 

I2/c2-Q2  =  2NkT  ,  (1.1) 

where  I  is  the  total  current  in  the  equilibrium,  Q  is  the  total  charge  per  unit 
length,  2NkT  is  the  total  thermal  energy  per  unit  length  and  c  is  the  speed  of 
light.  This  condition  can  he  understood  as  a  qeneral izat.ion  of  the  Bennett 
pinch  condition  (Reference  M. 
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The  quantities  on  both  sides  of  Equation  (1.1)  can  be  shown  to  be 
Lorentz  invariants.  Thus,  a  radial  equilibrium  in  one  frame  can  be  shown  to  be 
a  radial  equilibrium  in  all  other  Lorentz  frames.  This  property  of  Lorentz 
invariance  was  used  to  create  plasmoid  beam  equilibria  in  one  frame,  a  frame  in 
which  both  electrons  and  ions  moved  toqether  at  nearly  the  speed  of  light,  by 
setting  up  a  charge  neutral  (per  unit  length)  z-pinch  in  a  frame  where  electrons 
and  ions  moved  slowly  in  opposite  directions  at  equal  and  opposite  velocities 
p0  where  pQ<<1 .  in  this  z-pinch  frame  Equation  (1.1)  becomes, 

1  2/C^  =  i'^/C*  -  Q*  =  2N' kT ’  =  2NkT ,  (1.2) 


where  the  primed  quantities  are  transformed  into  the  z-pinch  frame.  The 
plasmoid  equilibrium  was  then  merely  a  rapidly  moving  z-pinch.  The  electron  and 
i o. i  densities  for  equilibrium  are  then  related  by  the  ratio 


n.  /n  =  y. h 
l  e  l  e 


(1.3) 


where  n|  and  ng  are  the  electron  and  ion  densities  respectively  and  where 
y  ^  and  ye  are  the  relativistic  gamma  factors  for  the  ions  and  electrons 
respectively,  y  =  (1  -  p2)*1/2  where  p  =  v/c. 


One  type  of  equilibria  that  was  found  was  that  of  a  Bennett  profile 
equilibrium  with  electrons  and  ions  having  similar  profiles.  The  profiles  were 
solutions  to  the  set  of  coupled  equations: 


dn  n  ue  x  f 

~kTe  =  —  J  "e^e  '  1}  (Ve  '  1)]  d* 
dx  \  o  y 


(1.4a) 


dn.  n.Tie^  x  y. 

—  kT  =  — -  J  n.[(p2  -  1)  -  —  (P.PA  -  D]  dx 

dx  x  o  Y . 


(1.4b) 


where  kTe  and  kT  ^  are  the  transverse  temperatures  of  electrons  and  ions, 
x  =  r^ ,  and  where  Pe  and  p^  are  the  normalized  electron  and  ion  velocities, 
respectively.  The  solutions  are  found  to  be 
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n  =  n  /(I  -  (x/a)V 
a  ao 


where  naQ  is  the  density  on  axis  for  either  specie  and  where  the  Bennett 
radius,  a,  Ls  defined  as: 


a 


z 


2kT  / 
a  a 

7iez( y.  Y  ( 1  -  P  p.  ) 

x  e  ex 


1 )  n 

ao 
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it  was  shown 
as  cosh  (e;- 
respect ively  , 

YekTe/ni’  so 
it  assumed  in 

i , 

a 


further  that  the  quantity  YiYfc  (1*PiPe)  could  be  expressed 
efe)  where  and  efc  are  equal  to  tanh"lp|  and  tanh'1 

and  that  this  quantity  was  itself  invariant  between  frames,  as  was 
that  the  Bennett  radius  in  all  frames  was  equal  to  the  value 
the  z-pinch  frame: 


kT  1 

ue^n’ p,z 
o  o 


(1.7) 


where  n1  and  kT  ’  are  the  common  densities  and  transverse  temperatures  of  the 
00  -i 

electrons  and  ions  in  the  z-pinch  frame.  The  quantity  becomes  equal  to 

cosh  (£|  -  ee)  in  the  limit  that  pQ,  the  drift  velocity  of  the  electrons 
and  ions  in  the  z-pinch  frame,  is  small.  The  slow  drift  speed  in  the  z-pinch 
frame  xs  one  of  rhe  initial  simpiifyinq  assumptions  of  the  model. 


It  has  been  observed  in  simulations  of  plasmoid  equilibria  of  general 
initial  profiles  that  electron  and  ion  profiles  evolve  very  quickly  to  the  same 
shape  if  their  initial  profiles  are  dissimilar.  Mso  it  has  bec-n  observed  that 
once  the  profiles  of  electrons  and  ions  become  similar,  there  is  a  slower 
evolution  to  a  Bennett  profile.  This  is  true  regardless  of  the  initial  profiles 
for  the  plasmoids. 


It  will  be  one  purpose  of  this  report  to  explain  why  similar  Bennett 
profiles  for  electrons  and  ions  should  be  expected  for  any  nonrotating  plasmoid 
after  many  ion  betatron  periods.  The  reasons  for  this  expectation,  it  will  be 
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shown,  involve  the  tendency  for  systems  to  seek  their  lowest  accessible  energy 
state  and  also  their  tendency  to  move  to  a  state  of  maximum  entrbpy .  Possible 
final  states  for  beams  with  initial  cannonical  angular  momentum  are  also  of 
interest  and  will  be  discussed. 

The  method  chosen  for  investigating  these  equilibria  issues  will  be  a 
variational  approach.  Variational  principals  have  been  shown  to  be  powerful 
techniques  for  predicting  the  end  conditions  of  very  complicated  systems.  They 
have  found  increasing  use  in  maqnet ic  confinement  fusion  theory  to  predict  the 
results  of  experiments  (Reference  IT. 

Because  of  its  compactness  we  will  use  fluid  theory  to  descrioe  beam 
equilibria  and  their  evolution.  The  same  results  could  be  obtained  through 
Vlasov  theory,  but  not  without  the  loss  of  simplicity  and  easy  understanding. 

1.4  BASIC  ASSUMPTIONS  OF  RELAXATION  MODELS 

When  a  system  evolves  to  maximize  or  minimize  some  global  quantity  the 
system  is  said  to  "relax".  The  relaxation  of  systems  to  minimum  energy  states 
or  maximum  entropy  is  observed  often  in  everyday  life.  The  movement  of  water  in 
a  sink  to  assume  a  common  lowest  level  when  more  water  Is  added  is  an  example  of 
enerqy  minimization .  The  diffusion  of  heat  in  a  frying  pan  until  its  surface  is 
a  uniform  temperature  can  be  thought  of  as  a  maximization  of  entropy.  It  seems 
reasonable  to  expect  similar  behavior  in  non-everyday  systems;  however,  in 
everyday  systems  the  mechanisms  for  relaxation  and  dissipation  are  clearly 
understood.  The  production  and  subsequent  phase-mixing  and  damping  of  water 
waves  allows  water  to  reach  a  common  level  and  the  process  of  heat  conduction  by 
microscopic  processes  causes  the  frying  pan  to  relax  to  a  uniform  temperature. 
Therefore,  when  we  consider  the  relaxation  of  intense  particle  beam  equilibria 
we  must  consider  what  process  or  processes  will  lead  to  dissipation  and 
increases  in  entropy.  In  beam  equilibria  the  particles  will  move  in  a 
complicated  system  of  interweaving  orbits.  Actual  close  approaches  by 
individual  particles  to  one  another  so  as  to  strongly  perturb  each  others  motion 
will  be  rare;  however,  stronq  oscillations  of  local  space-charge  density,  caused 
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by  th  collective  interact  ion  of  particles,  can  strongly  perturb  the  motion  of 
particles,  both  in  groups  of  part icles  passing  through  each  other  (a  t wo  stream 
instability )  at  hiqh  relative  speed  or  by  stronq  charqe  separations.  These 
oscillations  can  lead  to  rapid  exchange  of  energy  and  mutual  perturbations  of 
orbits  among  particles. 


This  in  turn  will  lead  to  a  more  uniform  spread  of  energy  among 

particles  and  a  more  diffuse  and  more  thoroughly  mixed  beam  profile.  The 
diffusion  amonq  particles  of  enemy  and  the  diffusion  in  space  of  particle 
orbits  by  electrostatic  and  electromagnetic  instabilities  can  cause  changes  in 

the  beam  profile  to  some  preferred  or  more  stable  state.  It  is  these  observed 
inst abil it ies  that  will  be  assumed  to  lead  to  relaxation.  These  oscillations 

can  be  expected  to  have  frequencies  of  the  order  of  the  beam  plasma  frequencies 

or  the  particle  betatron  frequencies.  The  frequencies  are  characteristic  of 
collective  particle  motion  in  the  fields  of  the  beam 


2elp, 

V  M  a*"c 
a  a 


(1.9) 


where  Ma  is  the  mass  of  the  particles. 


Oetailed  analysis  of  the  instabilities  is  very  complicated  for  general 
beam  profiles  and  will  not.  be  necessary  for  our  purposes.  We  will  merely  assume 
that  they  will  exist,  will  create  disorder,  and  dissipate  energy.  it  will  be 
assumed  that  the  oscillations  form  waves  that,  will  move  outward  and  be  lost  at  a 
very  large  radius.  The  loss  of  information  and  enerqy  from  the  system  will  be 
considered  irreversible. 


Because  we  are  concerned  with  the  evolution  of  radial  profiles  we  will 
assume  the  beam  equilibria  are  initially  uniform  in  the  z-direction  and  that 
this  condition  is  maintained  throuohout.  the  equilibrium  evolution  to  a  final 
state.  We  will  thus  not  allow  wave  propagation  in  the  z-direction.  This  wi 
have  the  effect  of  not  allowing  two  stream  instabilities  to  grow  from  the 
count  erst  reaming  of  the  electrons  and  ions  in  the  z-airection.  Because  this 
oounterstreaminq  is  assumed  to  be  weak  and  thus,  not.  a  large  source  of  free 
enerqy,  this  should  not  be  a  had  approximat ion . 


Certain  quantities  are  always  conserved,  even  in  complicated  systems, 
when  a  system  relaves.  Such  quantities  as  total  charge  per  unit  length,  current 
in  the  z-direction  and  particle  number  (per  unit  length)  will  be  assumed  to  be 
unchanged  b>  the  relaxation.  The  conservation  of  current  being  consistent  with 
our  approximat ion  of  uniformity  and  non-propaqat ion  of  waves  in  the  z-direction. 

Therefore,  we  will  assume  that  collective  oscillations  on  the  order  of 
the  betatron  frequencies  of  the  particles  will  provide  a  mechanism  for 
increasinq  entropy  in  the  beam  and  removinq  fee  enerqy  .  The  beam  will  thus 
minimize  its  free  enerqy  and  maximize  its  entropy .  However,  the  final  state  of 
the  beam  will  have  the  same  maqnet ic  enerqv  per  unit  lenqth. 

1.5  OVERVIEW 

In  the  second  chapter  of  this  report  the  theory  of  relaxation  of 
plasmoid  equilibrium  in  two-dimensions  will  be  discussed.  It  will  be  shown  that 
this  relaxation  can  be  studied  in  ail  frames  through  Lorentz  invariant  energy 
constraints,  that,  is  equivalent  to  a  relaxation  in  the  charge  neutral  z-pinch 
frame.  Because  of  this  the  z-pinch  frame  will  be  used  for  the  majority  of 
theory  and  simulations.  In  the  third  chapter,  results  of  simulations  will  be 
compared  with  theory.  A  final  chapter  will  summarize  the  results  of  this 
research . 
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SECTION  2 

THE  THEORY  OF  RELAXED  PLASMOID  EQUILIBRIA 


2.1  RELAXATION  IN  THE  Z -PINCH  FRAHt 

The  simplest  place  to  demonstrate  relaxation  of  a  plasmoid  to 

equilibrium  is  in  the  z-pinch  frame.  However,  the  basic  model  of  plasmoid 
equilibrium  requires  Lorentz  invariance,  so  that,  a  complete  theory  of  relaxat  ion 
should  apply  in  all  frames. 

One  of  the  key  assumptions  that  one  makes  in  the  z-pinch  frame  to 
obtain  an  equilibrium  is  the  assumption  of  similar  density  profiles  for 

electrons  and  ions.  In  the  z-pinch  frame  this  assumption  is  easy  to  justify  on 
the  grounds  that  this  minimizes  the  electric  field  energy,  since  all  charge 
densities  vanish  in  the  z-pinch  frame  for  similar  profiles.  However,  the 
assumption  of  similar  profiles  does  not  result  in  the  vanishing  of  electric 
fields  in  other  Lorentz  frames  because  the  charge  densities  are  shifted  by 

differing  gamma  factors.  Therefore,  one  way  to  make  the  Lagrange  minimization 
Lorentz  invariant  is  to  find  out  if  similar  profiles  for  electrons  and  ions  are 
required  for  an  equilibrium  or  are  just  a  convenient  assumption.  To  do  this  we 
will  reexamine  the  variational  process  in  the  z-pinch  frame  for  a  non-rotating 
plasmoid  and  see  if  anything  crucial  has  been  left  out. 

The  variational  problem  that  applied  in  the  z-pinch  frame,  where 

electrons  and  ions  were  assumed  to  have  similar  profiles,  could  be  written: 

D 

F  =  J  - —  +  \nln(n)  2mrdr  ,  (2.1) 

where  B^/Sn  is  the  magnetic  energy  and  n  is  the  density  and  nln(n)  can  be 
considered  a  thermody namic  entropy.  The  quantity  1  is  a  Lagrange  multiplier 
(Reference  2).  The  condition  of  maqnetic  energy  being  minimized  while  entropy 
is  kept  constant  yields  the  condition, 


CD 


OF  -  oAxB|  +  jo34  +  on  xin(n)  2ardr  ,  (2.2) 

I  z 

o 

*here  Az  is  the  vector  potential  and  where  we  consider  =  eon  v^  and  we 
consider  vz  constant . 

This  yields  the  equation 

ev4+\ln(n)=0,  (2.3) 

Z  Z 

and  its  derivative 

ev  nB  =  kT  —  .  (2.4) 

z  o  , 

dr 

2.2  THE  LOREN TZ  GENERALIZED  VARIATIONAL  PRINCIPLE 

This  last  equation  is  mereLy  the  equation  for  a  Bennett  pinch 

(Reference  3)  without  rotation.  The  generalized  variational  principle  takes  the 
form 

E"-8" 

F  -  -  +  a.  niK  i  in(n.kT.)  +  az  n  kT  in  (n  kT  )  (2.5) 
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Su 

+  kun?  -  (n.P.)^J  *  a4  ^n^  -  (n  p  ) 2 J  +  a  (p^  -  d2/c2)  , 
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where  E  is  the  electric  field,  Y|  and  ye  are  the  electron  and  ion  gamma 

factors,  respectively,  p[  and  pe  are  the  electron  and  ion  velocity 

normalized  to  c,  and  p2  -  d^/c2  is  the  electromagnetic  four  current.  This 

generalization  means  that  the  Lorentz  invariant  form  of  electromagnetic  field 
strength  is  minimized  while  Lorentz  invariant  entropies  and  densities  for  both 
species  and  the  total  Lorentz  invariant  charge  density  squared  are  all  kept 
constant.  The  results  of  this  variation  of  both  ne  and  n^  independently, 
with  the  usual  vanishing  of  surface  terms,  can  be  expressed  as: 
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4  2e  (n.-n  -n.p.^-n  p  p.)  =  0 
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(2.6) 


-e(  •?  -  A  p  )  +  a.  kT  In  n  kT  +  a.  2(n  -2pzn  ) 
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a...  2e  (n  -n,  -n  p  z-n.p  p. )  =  0 
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(2.7) 


The  presence  of  cross  terms  in  each  equation  requires  ny  ye  = 
ney :  and  thus,  similar  profiles.  Similar  profiles  for  electrons  and  ions, 
thus,  are  required  bv  minimization  of  L^-Bz  and  simultaneous  conservation  of 
particle  number  and  total  charge.  Differentiation  leads  to  the  equilibrium 
equation  Equations  1.4a  and  1.4b. 


ROTATIONAL  EQUILIBRIA  AM)  THE  GALACTIC  MODE 


The  inclusion  of  rotation  in  the  formalism  is  accomplished  by  the 
rotation  potential  for  each  specie 


£>.  =  m .  j 
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w  =  m  j 
e  e  1 
o 


dr' 


(2.9) 


where  and  mfc  and  VQ  -L  and  Vtte  are  ion  and  electron  masses  in  the  frame 

of  interest  and  rotational  velocities  respectively.  It  should  be  noted  that 
ni?i  and  ne-®e  dre  Lorentz  invariant.  Inclusion  of  strict  angular 

momentum  conservation,  it  should  be  noted,  results  in  unphvsical  equilibrium 
when  coupled  with  nonzero  temperature;  only  uniform  density,  cosmos  filling 
plasmoids  are  allowed.  However,  inclusion  of  nj^  and  ne<±>e  will  result 
in  a  general  class  of  rotating  equilibria. 


The  galactic  rotation  mode  results  when  the  condition 

4  x  .  a  i 
z  l  e 


(2.10) 
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is  applied.  This  condition  occurs  when  the  equilibrium  is  required  to  be 
insensitive  to  chanqes  in  the  entropy  or  temperature  of  the  species.  The 
rotational  velocity  profile  which  results  in  Fiqure  2  is  similar  to  that  seen  in 
qalaxies  (Reference  5).  Some  velocity  profiles  seen  in  galaxies  are  shown  in 
Fiqure  3.  The  qalact ic  mode  velocit>  profile  seen  in  Fiqure  2  can  be  written 
as, 


V  (r/a) 

V  =  — -  ,  (2.11) 

o  -  ’ 

v ,  c  ,  l 
1+r  /a 


where  V0  is  an  asymptotic  equilibrium  velocity.  The  conditions  for  the 
derivation  of  the  galactic  mode,  an  insensitivity  to  temperature,  indicate  a 
condition  of  nearly  circular  particle  orbits.  For  the  more  realistic  case  of 
very  noncircular  orbits,  the  conservation  of  angular  momentum  for  individual 
particles  will  cause  a  sharp  peak  to  occur  in  the  profile  at  smaller  radii 
as  particles  speed  up.  This  can  be  included  by  an  additional  term  in  Equation 
2.11. 

V  (1+r/a) 

ve  =  °~ .  r>a  .  (3.12) 

/1+r 2/a2 

V  =  /7  V  r/a  r<a  .  (2.13) 

6  o 

2.4  AXIALLY  MAGNETIZED  EQUILIBRIA 

The  case  of  rotational  motion  that,  produces  current,  and  from  this  axial 
magnetization,  can  be  examined  by  including  a  Lorentz  invariant, 

J  V1  -  2urdr  (2.14) 

while  eliminating  the  invariants  associated  with  entropy.  It  has  been  seen  that 
the  introduction  of  axial  magnetic  flux,  of  comparable  energy  to  thermal  energy, 
tends  to  cause  a  relaxation  that  is  determined  by  field  quantities  rather  than 
matter:  for  this  reason  thermodynamic  entropy  seems  to  be  unimportant. 
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StCTION  3 

SIMULATIONS  OF  STLADY-STATL  AND  RELAXING  EQUILIBRIA 


3.1  BASIC  REQUIREMENTS  FOR  SIMULATIONS 

Simulations  of  the  three  types  of  equilibria  discussed  were 
performed.  The  simulations  consisted  of  equilibria  initialized  in  steady-state 
and  nonequ '  1  ibriu.m  conf  iqurat  ions  that  were  allowed  to  relax  to  equilibrium. 

The  initializinq  of  steady-state  fields  and  currents  in  two-dimensions 
in  an  electromagnet ical ly  self-consistent  manner  proved  more  difficult  than 
anticipated.  However,  once  numerical  integrations  over  the  actual  particle 
charge  densities  and  current  were  used,  instead  of  idealized  analytic  ones, 
simulations  could  be  started  without  severe  electromagnetic  turbulence.  Because 
of  the  fully  electromagnetic  nature  of  the  code  MAGIC  (Reference  7),  the 
limitation  on  the  timestep  was  determined  by  the  Courant  condition 

At  >  rr  £  ,  (3.D 

c 

where  l  is  the  dimension  of  the  smallest  cell.  This  required  thermal  speeds  for 
ions  and  electrons  tnat  were  close  to  the  speed  of  light 

V  =  c/3  .  (3.2) 

t  h 

This  allowed  significant  particle  dynamics  to  occur  within  a  reasonable  number 
of  timesteps. 

3.2  SIMULATIONS  OF  NONROTATING  EQUILIBRIA 

Simulations  were  performed  in  the  z-pinch  frame  for  reasons  of 
simplicity.  Simulations  of  nonrotating  equilibria  were  performed  first.  The 
Bennett  profile  was  observed  to  be  stable  for  many  thermal  transit  times: 


where  m  Ls  the  ion  or  electron  mass,  n0  is  the  central  number  density,  and  p 
is  the  drift  velocity  normalized  to  the  speed  of  light.  A  level  of  fluctuations 
of  about  10%  of  the  mean  magnetic  field  and  particle  density  was  observed  even 
in  steady -st ate ,  as  can  he  seen  in  Figure  4.  Electrostatic  oscillations  were 
observed  to  be  strongest  in  the  center  of  the  pinch. 


A  nonequiiibr ium,  nonrotating  beam  was  also  simulated  and  allowed  to 
relax  to  equilibrium.  The  equilibrium  appeared  to  approach  the  Bennett  state 
after  many  thermal  transit  times.  Density  profiles  of  the  initial  and  final 
states  are  shown  in  Fiqure  5.  This  period  of  diffusion  was  much  longer  than 
expected.  The  expected  time  for  significant  diffusion  to  occur  was  a  few 
thermal  transit,  times. 


3.3  ROTATING,  NON  AXIALLY  MAGNETIZED  EQULIBRIA 


The  case  of  rotating  equilibria  was  simulated  by  being  initialized  in 
both  steady-state  and  nonequilibrium  configurations.  A  steady-state  rotational 
velocity  profile,  consistent  with  a  Bennett-like  profile,  was  arrived  at  by 
including  the  effect  of  large  radial  excursions  of  particles.  This  modified  the 
originallv  predicted  galactic  mode  as  is  shown  in  Figure  6.  The  rotating 
nonequilibrium  plasmoid  was  observed  to  relax  to  a  rotation  and  density  profile 
close  to  the  previously  simulated  steady-state. 

3.4  AXIALLY  MAGhCT IZED  EQUILIBRIA 

The  case  of  rotation  with  axial  magnet  izat,  ion  was  simulated  in  both 
steady-state  and  nonsteady -state  cases.  An  approximately  Bennett-  or 
Bessel -funct ion-shaped  profile  with  a  modified  galactic  mode  was  observed  for 
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Initial  (A)  and  final  (B)  maqnetic  field  profiles  for  a  Bennett 
plasmoid.  Not'-  larne  fluctuations  visible  in  final  state.  Sue 
fluctuations  occurred  throughout  the  simulation  and  appeared  to 
be  a  feature  of  the  equilibrium. 


.0 


to 


0.0 

radius 


radiu*  in  meters 


radius  in  meters 


«-  velocits  versus  radius  from  MAGIC  simulations  at  (A)  initial 
time,  and  (B)  after  ten  particle  transit  times,  for  an 
equilibrium  plasmoid.  The  "galactic"  profile  has  been  preserved 
with  onl\  a  smaLl  amount,  of  diffusion. 


the  steads -state  rotation  profiles  in  both  cases.  Bs  shown  in  Figure  7,  the 
initial  o  and  axial  magnetic  field  profiles  for  an  approximate  stead) - st at e 
profile  were  preserved  after  man)  particle  transit  times  with  only  minor 
changes.  The  final  state  of  a  nonst ead) - st  ate  initial  profile  also  appeared  to 
relax  to  this  state,  though  the  initially  nonst eady -st at e  profile  was  obseryed 
to  have  regions  of  axial  magnet  ic  field  reversal  in  late  times.  in  general,  the 
axiall)  magnetized  state  took  longer  to  relax  than  those  states  without  axial 
magnet izaton . 

3.5  RLLAXAT ION  TUCS 

Diffusion  from  uniform  profiles  to  Bennett  or  Bessel  function  stead)-state 
profiles  took  much  Longer  than  the  predicted  thermal  transit  time.  instead,  the 
diffusion  timescale  appeared  consistent  with  a  Bohm  diffusion  time  (Reference 
5 )  ; 


op 


(3.5) 


This  is  on  the  order  of  TOO  thermal  transit  rimes.  To  test  this  a  simulation 
was  run  with  the  same  Bennett  radius,  but  with  higher  drift  velocity,  thus 
giving  a  shorter  Bohm  time.  Signs  of  enhanced  diffusion  were  observed:  however, 
much  longer  simulations  would  have  to  be  run  to  give  d  finitive  proof  that  the 
Bohm  time  in  fact  controls  the  diffusion  process. 


3-6 


w  ■■  ;  it :  r:  «  v  *  b  r  7  «  v  •*  v  '•  ■ 


r>.-»  V  > 


t  •" 


s.\  / 


Figure  7.  Initial  (A)  and  final  (B)  e  magnetic  field  profiles  for  an 

axially  magnetized  plasmoid  core  equilibrium.  Axial  magnetic 
field  profiles  for  initial  (C)  and  final  (D)  states  of  the  same 
equilibrium. 
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SUMMARY  AM)  FUTURL  WORK 
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The  proqram  of  theory  and  simulations  has  made  great,  proqress. 
Steady-state  equilibria  of  three  basic  types  have  been  identified  theoretically 
and  verified  with  simulations.  The  theory  that  a  Lorentz  invariant  variational 
problem  could  predict  relaved,  robust  plasmoid  core  equilibria  has  been 
verified.  In  addition,  the  relaxation  process  has  been  observed  to  occur  for 
these  three  cases  and  its  timescale  has  been  approximately  characterized .  Ihis 
proaress  can  serve  as  a  foundation  for  an  expanded  effort  to  define  plasmoids  in 
three  dimensions  as  well  as  address  the  problem  of  propagation  of  plasmoids 
across  magnetic  fields. 
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